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'"^^This  report  presents  a heuristic  for  solving  multi-echelon  multi- indentured 
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that  consume  considerable  amounts  of  computation  time  and  the  procedures  are 
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determine  if  the  heuristic  can  also  be  extended  to  more  echelons,  lo 
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1.  INTRODUCTION 


In  1973,  Muckstadt  [2]  extended  Sherbrooke's  well  known  METRIC 
model  [3]  to  explicitly  account  for  multi- Indentured  Items.  By  a 
reparable  multi-indentured  Item  we  mean  a reparable  module  that  contains 
reparable  components.  The  components  themselves  may  contain  reparable 
subcomponents,  etc.  Much  of  the  Army's  complex  new  equipment  is  designed 
with  these  indenture  levels.  The  basic  idea  behind  indentured  items  is 
that  when  a module  fails,  a failed  component  may  be  quickly  removed  and 
replaced  with  a serviceable  component.  Thus,  the  actual  downtime  of  the 
module,  that  is,  the  time  the  module  is  not  in  a serviceable  condition,  may 
be  less  than  if  repair  had  to  be  done  on  the  whole  module.  Muckstadt 's 
MODMETRIC  model  leads  to  a complex  solution  technique  that  consumes  a con- 
siderable amount  of  computer  time.  In  this  report  we  present  a heuristic  for 
reducing  the  time  to  solve  MODMETRIC  while  yielding  close  to  optimal  solutions 

Section  2 presents  a brief  review  of  the  MODMETRIC  formulation  while 
in  Section  3 we  discuss  the  solution  procedure  in  more  depth.  Sections 
4 and  5 discuss  the  heuristic  we  propose  and  other  possible  heuristics. 

In  this  report  we  consider  a two  echelon  system  as  in  Figure  1. 


FIGURE  1 

There  are  N bases  and  a depot.  The  dynamics  of  a multi-echelon  multi- 
indentured  system  are  depicted  in  Figure  2.  We  will  consider  a module  with 
M components.  When  referring  to  the  module  and  its  components  a 0 subscript 

will  refer  to  the  module;  the  components  will  be  numbered  1,  2 M. 
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- FAILED  MODULE 
"FAILED  COMPONENT 
• SERVICEABLE  MODULE 
^•SERVICEABLE  COMPONENT 


2.  THE  MODMETRIC  PROBLEM 

In  this  section  we  present  a brief  review  of  the  MODMETRIC  problem 
as  formulated  by  Muckstadt  [2]. 

The  basic  assumptions  are: 

(1)  Customer  demand  at  each  base  for  modules  is  a stationary 
(compound)  Poisson  Process  with  rate  X^ . 

(2)  No  lateral  resupply  among  bases. 

(3)  Failures  of  one  type  of  item  (i.e.  modules  or  components) 
are  Independent  of  failures  of  other  types  of  items. 

(4)  All  repair  times  for  all  items  are  statistically  Independent. 

(5)  The  echelon  at  which  an  item  is  repaired  depends  only  on  the 
complexity  of  the  repair  required. 

(6)  No  batching  (or  waiting)  of  items  before  repair  is  begun. 

(7)  All  locations  follow  a one-for-one  (S-1,S)  resupply  policy. 

For  ease  of  exposition  we  further  assume: 

(8)  All  parts  are  repaired  (no  condemnations) . 

(9)  A module  failure  is  caused  by  the  failure  of  at  most  one 
component . 

(10)  Modules  sent  to  the  depot  for  repair  are  completely  overhauled. 

Assumptions  (8)  and  (10)  are  not  restrictive  and  could  easily  be 
relaxed  in  the  model  formulation.  In  actual  Army  depot  situations  a 
module  is  sent  to  the  depot  by  a base  only  when  the  module  has  been  seriously 
damaged.  In  these  situations , the  time  to  overhaul  the  module  is  large 
enough  so  that  any  component  that  is  removed  from  a module  is  repaired 
and  the  now  aai’&tceable  component  is  put  back  into  the  same  module  from 


w gaffi 

Sa  remc 


which  it  w Se  removed. 

f 

The  problem  we  consider  is  to  determine  the  module  and  component 
stock  levels  at  all  locations  that  minimize  the  sum  of  the  holding  cost 
of  all  items  (modules  and  components)  at  all  locations  (bases  and  depot) 
subject  to  a performance  constraint  on  the  sum  over  all  bases  of  the  time 
weighted  module  backorders. 

In  order  to  compute  the  expected  number  of  module  backorders  out- 
standing at  a point  in  time  at  base^ , we  need  to  know  the  module  resupply 
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time  at  baae^ . It  is  through  the  expected  module  resupply  time  that 
MODMETRIC  explicitly  accounts  for  the  effect  of  the  components  on  module 
performance. 


Let  Tj  - expected  module  resupply  time  at  base^ • 


Then 

(1) 

Where 


‘s 


VbJ  + V + (1  ' V (°3  + V 


proportion  of  module  failures  repaired  at  base 


j 


Let  S 


b * mean  base,  fault  Isolation  and  component  remove  and 
J replace  time 

Z.  ■ expected  delay  In  module  repair  at  base,  due  to 
J unavailable  serviceable  components  ^ 

Oj  - mean  order  and  ship  time  between  base^  and  depot 

expected  delay  at  depot  due  to  unavailable  modules 

spare  stock  level  of  item  k at  location  j;  k-0, 
1,...,M  (0  ■ module);  j ■ 0,1,..., N (0  s*  depot) 


'j 

kj 

*3 


expected  dally  module  demand  at  base 


j 


Xq  ■ expected  daily  module  demand  at  depot 


N 

- I 


3-1 


yi-v 


F.  ■ mean  depot  repair  time  for  item  k;  k 
1, . . . ,M  (0  - module) 


0, 


The  expected  delay  at  the  depot  from  the  time  base^  places  a module 
resupply  request  on  the  depot  until  a serviceable  module  is  available  for 
shipment  to  base^  is  given  by 


expected  number  of  module  backorders  out- 

D.  >*  standing  at  a point  In  time  at  depot 

J expected  dally  demand  for  modules  at  depot 


Given  a depot  module  stock  level  of  Sqo  the  expected  module  back- 
orders outstay  ^ng  at  the  depot  at  a point  In  time  la  given  by 


(2) 

where 

Hence , 
(3) 


x > S 


(x  - Soo)  p(x;  XJJ 


oo 


o o 


p(x;u. ) - probability  that  x units  are  In  resupply  at  location 
J j given  a mean  resupply  time. 


x > S 


(x  - S)  p(x;  X F ) 
00  o o 


oo 


Dj  reflects  the  Interaction  of  the  two  echelon  supply  systen  on  module 
performance  at  base^ . Similarly  Z^  reflects  the  effect  of  the  components 

on  module  performance  at  base^ . We  now  show  how  may  be  calculated. 

Let  P.  - Probability  that  a module  failure  was  caused  by  a failure 
of  component  k k“l M. 

g.  - Expected  delay  In  base,  module  repair  given  that  the 

J module  failed  due  to  aJ failure  of  component  k. 


kj 


Dally  demand  for  component  k at  location  j 
k*l, ... jM  j*0> ... jN 


Then,  clearly 


(4) 


M 

Z . - Z 


Since  module  failures  that  require  base^  repair  follow  a Poisson 
Process  with  rate  r^X^  then  component  k failures  follow  a Poisson  Process 
with  rate  X. 


kj 


Hence, 

(5) 


Vjrj 


X, 


P - ^ 
k Vj 


l 


To  calculate  we  need  to  know  the  base^  resupply  time  for  component 
If  T.  , - base . resupply  time  for  component  k,  then  similar  to  (1) . 


where 


Proportion  of  component  k failures  repaired  at  base 


Mean  component  k repair  time  at  base 


Component  k order  and  ship  time  between  depot  and  base 


Delay  at  depot  due  to  unavailable  serviceable  component 
k stock. 


Given  a depot  spare  stock  level  of  of  component  k,  we  have,  as 
before,  that 

. E , <x  - Sko>  O'*  \„V 

X > 0| 


Substituting  (5)  and  (8)  Into  (A)  we  have  finally  that 


With  this,  all  the  terms  In  the  expected  base^  module  resupply  time 
equation  (1)  have  been  calculated. 


In  summary  then,  given  the  system  parameters  and  module  and  component 
stock  levels  at  all  locations,  to  calculate  the  base^  module  resupply 
time,  we: 

(a)  calculate  the  delay  at  depot  due  to  unavailable  components 
of  type  k;  k - 1,...,M. 

(b)  calculate  the  component  resupply  time  at  base,  for  each 

component . ^ 

(c)  calculate  expected  component  backorders  at  base  and  then 

the  delay  In  base^  module  repair  due  to  unavailable  components. 

(d)  determine  expected  delay  at  depot  due  to  unavailable  modules. 

(e)  calculate  base  module  resupply  time. 

It  should  be  clear  from  the  above  discussion  that  to  explicitly 
model  the  component  - module  relationship  leads  to  some  rather  complex 
expressions  and  a difficult  optimization  problem.  In  the  next  section 
we  explore  solution  procedures  for  the  optimization  problem  in  more  depth. 


8 


3.  OPTIMIZATION  PROBLEM 

The  MODMETRIC  Optimization  Problem  can  now  be  formulated  as 


N M 

Min  E C°  S + E cj  S 

j-0  H °j  k-1  H kJ 


N 

S.T.  E 


j-1  x > Sr 


(x  - S0j)  p(x;  \i  Tj)  £ 


Where  £ ■ (s00*  Soi’***’  Son:  sm sin»  s->n smn^ 


IN’  20' 


is  a vector  of  non-negative  integers,  y„  is  a specified  module  performance 

k M 

level  and  C„  is  the  unit  price  for  item  k,  k-0,...,M. 
n 

Problem  (10)  is  not  convex  and  optimization  is  difficult  due  to  the 
complex  interactions  between  echelons  and  between  modules  and  components  as 
expressed  in  the  resupply  time  equations.  Kotkin  [1]  has  developed  bounds  on 
the  optimal  module  and  component  stock  levels  (and  hence  bounds  on  the  optimal 
module  and  component  budget  expenditures).  Shay  and  O'Malley  [4]  have  developed 
a solution  procedure  that  is  not  as  complex  as  the  procedure  suggested  by 
Muckstadt.  Even  so,  no  relatively  easy  solution  to  (10)  is  known. 

Muckstadt  [2]  suggests  partitioning  (10)  into  2 subproblems: 


(11a) 


N M 

Min  E E S 

j-0  k-1  H 

S'T‘  jE-l  k-1  x l S <X  ’ Skj)  P(X;XkjTkj)  -YC 


Then,  after  obtaining  the  solution  to  (11a)  we  solve: 


(lib) 


hi-  r_o  qs0j 


N 

S.T.  E 


j-1  x > S. 


(x  - S ) p(x;X  T ) < yM 
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r 


1 


where  f is  the  module  resupply  time  at  base^  given  the  component  stock  levels 
determined  in  (11a) . Note  the  variables  in  (11a)  are  the  component  stock 
levels  while  in  (lib)  the  variables  are  the  module  stock  levels.  Note  that 
(11a)  and  (lib)  are  problems  of  the  METRIC  [3]  type.  Muckstadt  solved  the  dual 
problems  to  (11a)  and  (lib)  (the  dual  problems  being  to  minimize  the  sum 
over  all  bases  of  the  expected  item  backorders  outstanding  at  a point  in 
time  subject  to  a budget  constraint)  and  suggested  the  following  procedure 
for  solving  problem  (10) . Set  lower  and  upper  bounds  on  component  invest- 
ment and  a component  budget  increment.  Then  for  a fixed  total  budget 
solve  the  dual  of  (11a)  with  the  component  budget  set  at  its  lower  bound 
and  then  solve  the  dual  of  (lib)  with  the  module  budget  - total  budget  - 
component  budget.  Next,  Increment  the  component  budget  and  repeat  the 
procedure.  The  approximate  solution  to  (10)  is  then  the  best  mix  of 
component  and  module  budgets  determined  by  this  procedure. 

The  procedure  Muckstadt  suggests  for  solving  (10)  can  involve  solving 
many  subproblems  of  the  type  (11a)  and  (lib) . This  procedure  may  consume 
a considerable  amount  of  computer  time  and  leads  to  long,  complex  computer 
programs . 

The  purpose  of  this  report  is  to  present  a heuristic  procedure  for 
solving  (10)  that  Involves  solving  subproblems  (11a)  and  (lib)  only  once 
each.  This  savings  in  solution  complexity  and  running  time  is  especially  de- 
sirable since  in  most  practical  situations  a tradeoff  curve  of  budget  versus 
performance  is  desired  rather  than  a solution  to  (10)  for  a particular  per- 
formance level  yu.  To  construct  this  trade  off  curve  involves  solving  (10) 
with  various  values  for  and  thus  the  need  for  a good  heuristic  for  solving 
(10)  is  apparent. 

By  introducing  a generalized  Lagrange  multiplier  in  (11a)  and 
multiplier  in  (lib)  these  problems  can  be  rewritten  as: 


(12a) 


Min 


N 

£ 

J-l 


M 

£ 

k-1 


<CH  Skj  + BC 


M 

£ 

k-1 


£ 

> 


3kj 


(x  ■ V p(x;XkjTkj)3 


c£  sv 
H ko 


and 
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— , 

ICS  s0j  + “m  J <»  - V V1  + <£  Soo 

x Oj 

Note  (12a)  Is  separable  by  component.  Furthermore,  the  multiplier  In 

(12b)  is  the  same  as  the  multiplier  we  might  use  in  (10)  if  we  were  to 

try  and  solve  (10)  directly.  By  specifying  B^  in  (10)  (and  hence  in  (12b)) 

we  are  implicitly  specifying  a target  module  performance  value  We 

then  ask  if  given  the  multiplier  Bu  can  we  a priori  determine  the  "optimal" 

* ^ 

multiplier  value  B^  in  (12a)?  If  this  can  be  done  then  solving  problems 
(12a)  and  (12b)  once  each  will  yield  a satisfactory  solution  to  (10) . The 
heuristic  we  present  will  be  an  attempt  to  find  the  "optimal"  B^,  given 

The  entire  approximation  procedure  that  we  employ  will  be: 

a.  Set  Module  backorder  penalty,  B^,  in  (10) . 

b.  Set  Component  backorder  penalty,  B^ . 

c.  Use  in  (12a)  to  determine  component  stock  levels. 

d.  Adjust  module  resupply  time  to  reflect  component  delays. 

e.  Use  B^  in  (12b)  with  adjusted  module  resupply  time  to  determine 
optimal  module  stock  levels. 

f.  Return  to  step  1 and  adjust  to  achieve  the  desired  target 
module  performance. 

This  paper  focuses  on  a heuristic  for  determining  B^  in  step  b.  above. 


(12b) 


N 

Min  I 

j-1 


r 


4.  HEURISTICS 

* 

In  che  next  section  we  present  s heuristic  for  determining  . 

Before  proceeding  with  this  heuristic  we  will  exsmlne  two  other  possible 
* 

heuristics  for  B^  . 

An  apparent  heuristic  would  be  to  set  Bj,  ■ ^ i.e.  charge  the  system 
as  much  for  a component  backorder  as  we  do  for  a module  backorder.  This 
would  imply  that  a component  backorder  is  considered  as  bad  as  a module 
backorder.  The  constraint  in  (10)  is  a constraint  on  module  backorders. 

But  not  every  component  backorder  Is  causing  a module  backorder.  If  a 
component  is  backordered  this  means  a module  Is  waiting  to  be  repaired. 

It  does  not  necessarily  mean  that  the  waiting  module  is  backordered.  So 
a component  backorder  may  only  mean  repair  on  a module  is  being  delayed 
and  in  this  case  we  would  not  want  to  charge  the  same  penalty  as  we  do 
on  a module  backorder.  Hence,  It  would  seem  that  Is  an  upper  bound 
on  the  component  backorder  cost.  By  charging  too  high  a component 
backorder  penalty  in  (12a)  we  would  be  investing  more  money  in  components 
than  we  would  normally  do  in  an  optimal  solution  to  (10) , thereby  causing 
an  under  Investment  in  modules. 

Note  that  since  B^  ■ in  this  heuristic  (called  the  heuristic) , 
Increases  the  heuristic  will  Increase  B^  and  thus  Increase  the  com- 
ponent budget  as  well  as  the  module  budget  to  meet  higher  module  performance 
targets. 

Another  possible  heuristic  would  be  to  set  B^  - module  unit  price 
(MUP  ■ Cy) . Since  a component  backorder  means  a module  is  waiting  to  be 
repaired  and  hence  Is  not  in  a serviceable  condition.  It  may  be  reasonable 
to  charge  the  system  the  unit  price  of  a module  that  is  no  longer  available 
to  It.  This  approach  however,  also  does  not  differentiate  between  the 
possible  effects  of  component  backorders.  Contrary  to  the  B^  heuristic, 
the  module  unit  price  heuristic  (MUP  heuristic)  says  component  backorders 
should  only  be  charged  for  delaying  repair  of  modules  and  not  for  causing 
module  backorders.  It  seems  reasonable  (and  we  shall  show  this  later) 
that  the  MUP  should  be  a lower  bound  on  the  true  value  of  the  component 
backorder  penalty  cost.  Underestimating  B^  results  In  an  underinvestment 
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in  components  thereby  causing  an  overinvestment  in  modules. 

Note  that  in  the  MUP  heuristic,  since  the  module  unit  price  does  not 
change,  component  stock  levels  do  not  change  for  different  values  of 
Although  this  heuristic  is  extremely  easy  to  implement  (since  it  involves 
solving  (12a)  only  once  and  these  component  levels  are  the  stock  levels 
for  all  values  of  B^)  it  seems  reasonable  to  expect  that  the  component 
levels  should  vary  for  different  values  of  the  module  backorder  penalty 
cost. 

Figure  1 through  Figure  4 illustrates  how  these  heuristics  compare 
with  Mucks tad t' s solution  procedure.  In  these  graphs,  we  have  plotted 
expected  backorders  versus  total  budget. 

The  four  curves  compare  the  and  MU?  solution  heuristics  with 
Muckstadt's  MODMETRIC.  From  these  figures  we  can  make  the  following 
observations: 

a.  For  low  ready  reates  (ready  rate  is  the  percentage  of  time  no 
module  backorders  are  outstanding)  both  the  B^  and  MUP  heuristics  work 
well.  (The  numbers  next  to  the  MODMETRIC  points  indicate  the  ready  rate) . 

b.  As  the  ready  rate  Increases  the  B^  heuristic  begins  to  do 
better  than  the  MUP  heuristic.  (Figures  1 and  2) . 

c.  As  the  ready  rate  Increases  even  further,  the  MUP  heuristic  begins 
to  do  considerably  better  than  the  B^  heuristic.  However,  both  are  generally 
inferior  to  the  MODMETRIC  solution  obtained  by  Muckstadt's  procedure 
(Figures  3 and  4) . 

Since  in  most  practical  cases  we  are  Interested  in  ready  rates  in 
the  90  - 100X  range.  Figures  3 and  4 imply  that  for  higher  ready  rates  it 
is  better  to  overinvest  in  modules  than  to  underinvest  in  modules.  This 
corresponds  to  an  empirical  observation  made  by  Muckstadt  [2]  in  his  original. 
MODMETRIC  paper. 

To  gain  some  insights  into  the  observations  we  shall  explore  (10) , 

(12a)  and  (12b)  in  more  depth. 

Claim:  If  the  module  stock  is  constrained  to  be  zero  at  all  bases, 

* 

then  Bc  - B^  for  all  B^. 

Proof:  By  introducing  the  multiplier  in  (10)  we  get  the  new 
problem. 
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FIGURE  1 
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x > S 


Min  C?S  + 
H oo 


N 

Z 

j-1 


oj 


(x-SQj)  p(x;XjTj)] 


If  Che  bases  are  constrained  Co  have  zero  module  stock  Chen  (13) 
reduces  Co: 

S M k 0 

d4)  Min  Z Z [C*  S + B X.T  ] + C„  S 

j-1  k-1  H M J J H oo 

Using  (1),  (3)  and  (9)  this  reduces  Co: 

N M M 

(15)  Min  Z [Z  [C„S  + K Z S (x  - S ) p(x;X..T  )] 

j-1  k-1  H kJ  * k-1  x > S.  . kJ  kJ  kJ 

+ bm  * (»-s0)P(xa0F0)  + C°  soo 

x>S 

o 

But  (15)  Is  precisely  problems  (12a)  and  (12b)  with  the  multiplier  value 
equal  Co  and  the  claim  Is  established. 

Intuitively,  If  there  is  no  module  stock  at  the  bases,  every 
component  backorder  Is  directly  causing  a module  backorder  since  the 
module  waiting  for  repair  must  be  due  out  to  a customer.  Hence,  a com- 
ponent backorder  Is  as  bad  as  a module  backorder  and  It  seems  reasonable 
to  charge  the  same  backorder  cost,  to  component  backorders  as  we  do 
to  module  backorders. 

Note  that  if  B^  MUP  then  it  never  pays  to  have  any  module  stock 

since  It  costs  less  to  have  one  unit  year  of  module  backorders  than  to 

hold  the  module  for  one  year.  Hence,  to  force  base  module  stock  to  be 

zero,  set  B^  - MUP.  Unless  the  module  has  a low  demand  rate,  zero  module 

stock  levels  at  the  bases  will  result  in  a relatively  low  ready  rate. 

Hence,  for  low  ready  rates  we  would  expect  to  see  that  the  two  heuristics 

agree  with  each  other  and  with  MODMETRIC. 

The  above  observations  help  to  illustrate  why  MUP  Is  a practical  lower 
* 

bound  on  B.  . In  most  problems  of  interest  Bu  > MUP  since  we  desire  the 
system  to  have  spare  modules,  and  since  we  would  expect  the  component 
backorder  cost  to  Increase  as  the  module  backorder  cost  Increases,  MUP 
should  be  a lower  bound  on  B^*. 

As  B^  Is  Increased  above  MUP  and  we  start  to  put  module  stock  at  the 
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bases,  the  two  heuristics  begin  to  give  different  answers.  Initially 
the  overinvestment  In  components  by  the  B^  heuristic  is  less  damaging 
than  the  underinvestment  in  components  by  the  MUP  heuristic.  This  Is 
probably  due  to  the  fact  that  for  low  ready  rates,  achieving  a module 
performance  target  can  be  done  by  investing  in  either  modules  or  components. 
However,  for  higher  and  higher  ready  rate  targets  we  would  expect  modules 
to  play  a more  significant  factor  in  achieving  the  ready  rate  target. 

Hence,  the  HUP  heuristic  which  buys  more  modules  than  the  heuristic 
will  tend  to  do  better  for  higher  ready  rates. 

A procedure  for  solving  (10)  similar  to  the  Muckstadt  solution  pro- 
cedure would  be  to  specify  a multiplier  increment  d£.  Then,  for  a particular 
value  of  B^,  solve  problems  (12a)  and  (12b)  with  B^,  ■ MUP,  Bc  - MUP  + dc, 

■ MUP  + 2d£, Bc  - B^.  The  solution  would  then  be  the  best  of 

these  trial  solutions.  Here  again,  however,  many  solutions  of  (12a)  and 
(12b)  are  required  for  a particular  value  of  The  next  heuristic  we 

present  will  attempt  to  take  advantage  of  the  knowledge  gained  from  the 
MUP  and  B„  heuristics. 


5.  THE  FINAL  HEURISTIC 


r 


! 


Previously,  we  mentioned  that  specifying  In  (10)  and  (12b) 
Implicitly  specifies  a module  performance  target.  Particularly,  regardless 
of  the  component  levels,  specifying  establishes  a lower  bound 

(16)  A - 1 - 


MUP 


on  the  module  ready  rate  at  a base.  To  see  this,  note  that  If  the  component 
levels  are  fixed  at  all  locations  and  the  depot  module  stock  Is  fixed, 
then  (10)  (or  12b)  reduces  to 

(x  - SoJ)  p(x;XjTj)] 

oj 

where  f reflects  the  fixed  component  and  depot  module  levels.  (17)  Is 

separable  by  base  and  for  each  base  (17)  Is  convex  in  S . (This  Is  well 

oj 

known  and  a proof  can  be  found  In  Sherbrooke  [3]).  The  problem  for  base 

3 is 


(17) 


Min 


N 

£ 

J» 


[c„  s 


H oj 


"h 


£ 

> 


(18)  Min  VV  ' C2  SoJ  + “m  2 <«  - V 

x > S , 

oj 

The  optimal  base  j module  stock,  sQj*»  must  satisfy  the  optimality 
conditions : 

(a)  Qj  (SQj*)  - Qj  (SoJ*  - 1)  < 0 

(19)  (b)  Q(SqJ*  + 1)  - Q (Soj*)  > 0 
(c)  or  SQj*  - 0 

Using  these  conditions  we  see  that  SQj*  must  satisfy 

C° 

(20)  P (SoJ.  + liljTj)  < J-  . ^ < P (..j-.Xjfj) 
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1 1 


where 


p (s  ;v)  « z * 

J k-S 

oj 


p(fc'.p) 


it  it 

?_  (Sqj  + 1)  Is  the  probability  that  (Sqj  + 1)  or  more  modules  are  in 
resupply  to  base  j.  But  this  is  Just  the  probability  that  there  is  at 
least  one  module  backorder  at  base  j and  hence  the  ready  rate  ■ 1 - 
probability  of  one  or  more  module  backorders,  is  bounded  below  by  (16). 

Given  that  a component  is  backordered,  it  is  possible  that  a module 
backorder  exists  as  well.  If  this  is  the  case  it  seems  reasonable  to 
charge  B^{  for  the  component  backorder.  If  there  is  no  module  backordered 
then  we  charge  MUP  for  the  component  backorder.  However,  A in  (16)  gives 
an  approximation  to  the  percentage  of  time  that  there  will  be  no  module 
backorders  and  hence  1 - A is  the  proportion  of  time  there  are  module 
backorders.  Combining  the  above  observations,  the  heuristic  we  propose 
is  to  set 


(21) 


Bc*  - A (MUP)  + (1  - A)  (Bm) 

- (1  - §E£)  MOP  + (jj^)  bm 


In  summary,  the  heuristic  says  that  if  on  a day  on  which  there  is 
a component  backorder  there  are  also  modules  backordered,  charge  per 
component  backorder  day.  Otherwise  charge  MUP  per  component  backorder  day. 

A underestimates  the  true  ready  rate  since  A is  only  a target  and  it 
is  usually  exceeded  (see  equation  (20)).  However,  the  fact  that  there  is 
a component  backordered  affects  (lowers)  the  probability  that  there  are  no 
modules  backordered  at  the  same  time.  Hence,  we  use  A as  an  approximation 
to  the  probability  of  no  module  backorders. 

The  heuristic  (21)  has  several  desirable  properties.  From  (21)  we 
note  that: 
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(a)  MUP  < Bc*  < ^ 

(b)  Bc*  Increases  as  Increases 

(c)  As  goes  to  Infinity,  Bc*  goes  to  2 MUP 

Properties  (a)  and  (b)  are  desirable  In  light  of  the  discussion  of 
the  previous  section.  There  seems  to  be  no  Intuitive  significance  to  the 
value  2 MUP.  As  Bj^  increases,  Bc*  approaches  this  value  of  2 MUP.  Hence, 
for  large  values  of  B^,  Bc*  at  worst  underestimates  the  true  optimal 
component  backorder  cost.  As  mentioned  previously,  for  higher  ready  rates 
It  Is  better  to  underinvest  in  components  than  to  overinvest. 

Figures  5,  6 and  7 compare  the  heuristic  and  MODMETRIC  results.  These 
are  again  plots  of  expected  backorders  versus  total  budget.  Tables  1,  2 
and  3 give  the  corresponding  numerical  values  for  the  points  plotted  in 
Figures  5,  6 and  7 respectively.  From  the  graphs  and  tables  we  note 
that: 

(a)  the  heuristic  generates  some  but  not  all  of  the  MODMETRIC 
point 8. 

(b)  when  the  heuristic  generates  points  that  are  not  MODMETRIC 
points,  the  heuristic  points  lie  on  or  near  the  extended 
MODMETRIC  curve.  This  extended  curve  is  simply  a straight 
line  completion  of  the  MODMETRIC  points. 

(c)  most  of  the  later  points  generated  by  the  heuristic  are  ob- 
tained by  Increasing  the  module  stock  levels  and  not  the 
component  levels.  This  Is  because  all  Bc*  for  these  later 
points  are  nearly  equal  to  2 MUP. 

(d)  for  very  high  ready  rates,  the  heuristic  underinvests  In 
components  and  hence  overinvests  In  modules. 

When  the  heuristic  and  MODMETRIC  points  agree,  both  the  component  and 
module  stock  levels  at  all  locations  agree.  In  principle,  Muckstadt's 
MODMETRIC  solution  procedure  will  generate  a solution  for  any  budget  value. 
Many  of  the  heuristics  points  could  be  generated  by  the  MODMETRIC  algorithm 
if  in  the  search  algorithm  the  component  budget  increment  is  made  small 
enough.  In  general,  however,  the  heuristic  will  generate  fewer  points  than 
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MODMETRIC  does.  This  Is  because  as  B * approaches  2 MUP  the  component 

v 

levels  stay  the  same  and  only  module  levela  change.  On  the  other  hand, 
MODMETRIC  may  change  the  component  levels  as  well  as  the  module  levels 
and  hence  may  generate  more  points. 

Figures  8 and  9 are  plots  of  component  budget  verus  total  budget. 

From  these  plots,  we  see  that  MODMETRIC  may  "oscillate"  between  component 
levels.  That  Is,  as  the  total  budget  Increases,  the  allocated  component 
budget  may  oscillate.  Generally  we  would  expect  that  as  the  total  budget 
Increases,  the  component  budget  Increases.  Say  we  have  a total  budget  of 
$T  and  given  this  total  budget  we  have  an  optimal  component  budget  of 
$C  >_  MUP.  If  we  Increase  the  total  budget  by  (N-l)  MUP  where  N ■ number 
of  bases.  It  may  pay  to  decrease  component  Investment  by  MUP  so  that  one 
additional  module  can  be  put  at  each  base.  This  Is  one  of  the  situations 
In  which  total  budget  may  Increase  but  the  component  budget  may  decrease. 
Generally,  MODMETRIC  "oscillates"  before  levelling  off  at  a particular 
component  budget  since  after  the  total  budget  Is  Increased  enough,  there 
will  be  enough  money  for  more  modules  and  more  components. 

The  heuristic  obviously  cannot  oscillate.  Since  Increases  as 
Increases,  the  component  budget  can  never  decrease  as  the  total  budget 
Increases.  The  heuristic  tends  to  reach  the  component  budgets  where 
MODMETRIC  levels  off,  but  it  cannot  generate  the  points  where  MODMETRIC 
oscillates.  The  heuristic  may  generate  the  total  budgets  Tq  and  (Tq  + N MUP) 
while  MODMETRIC  may  generate  points  Tq,  (Tq  + (N-l)  MUP)  and  (Tq  + N MUP). 
This  helps  to  explain  why  MODMETRIC  will  usually  generate  more  points  than 
the  heuristic. 

The  figures  5,  6,  7,  8 and  9 are  typical  of  the  results  we  had. 
Experience  has  Indicated  that  the  heuristic  yields  points  satisfactorily 
close  to  the  MODMETRIC  points  to  be  a viable  alternative.  Note  that  the 
heuristic  (21)  adjusts  Itself  for  different  priced  modules  and  for  varying 
system  parameters. 

A good  decision  making  policy  might  be  as  follows.  Use  the  heuristic 
to  generate  the  performance  versus  budget  tradeoff  curve.  Then  a manager 
may  choose  the  point  he  wishes  to  be  at.  Once  this  point  is  determined. 
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the  more  precise  solution  procedure  (or  the  MODMETRIC  solution  procedure) 
may  be  used  to  obtain  a more  exact  answer  for  the  chosen  point.  However, 
even  for  the  chosen  point,  experience  indicates  the  heuristic  solution  will 
be  satisfactory  without  more  complex  and  time  consuming  procedures. 
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CONCLUSIONS 


In  this  paper  we  have  presented  a heuristic  to  solve  the  MODMETRIC 
inventory  problem.  Computational  experience  has  indicated  that  the 
heuristic  works  extremely  well.  The  heuristic  leads  to  a much  less 
complex  computer  code.  Furthermore,  the  heuristic  will  generate  a budget 
versus  performance  tradeoff  curve  at  least  ten  to  twenty  times  faster  than 
the  MODMETRIC  procedure.  In  Appendix  I we  have  a listing  of  the  computer 
code  for  the  heuristic. 

F-rther  research  will  be  done  to  determine  if  the  heuristic  can  be 
extended  to  more  than  two  echelons  and  more  than  one  level  of  Indenture. 
Preliminary  indications  seem  to  indicate  that  the  heuristic  is  easily  and 
successfully  extendable  to  more  levels  of  indenture.  The  heuristic  for 
subassemblies  is  just  modified  to  reflect  the  unit  price  and  the  backorder 
penalty  of  the  next  higher  assembly. 

Work  is  now  underway  to  extend  the  heuristic  for  problems  with  more 
than  two  echelons. 
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APPENDIX  II 


The  test  module  used  to  verify  the  heuristic  had  the  following 
characteristics : 

Module  unit  price:  $80,000 

Mean  hours  between  module  demands:  260 

Depot  module  repair  time:  60  days 


All  components  required  depot  repair.  The  component  data  was: 


Component  if 

Mean  Hours 

Unit 

Depot 

Between 

Price 

Repair 

. 

Demands 

Time 

1 

1000 

25000 

45 

800 

1000 

45 

3 

3800 

35000 

45 

4 

4000 

1100 

45 

5 

2400 

30000 

45 

6 

2200 

1500 

45 

The  order  and  ship  time  between  bases  and  depot  was  15  days  for  all  Items. 
There  were  two  identical  bases.  We  assumed  4,  8 and  12  module  failures 
per  base  per  month  in  Figures  5,  6 and  7 respectively. 

Other  cases  with  different  unit  prices,  repair  times,  failure  rates, 
etc.,  were  run  and  the  results  were  comparable  to  those  presented  in  the 
report . 
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